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We investigate a link between the running of the fine structure constant a and a time evolving 
scalar dark energy field. Employing a versatile parameterization for the equation of state, we 
exhaustively cover the space of dark energy models. Under the assumption that the change in a is to 
first order given by the evolution of the Quintessence field, we show that current Oklo, Quasi Stellar 
Objects and Equivalence Principle observations restrict the model parameters considerably stronger 
than observations of the Cosmic Microwave Background, Large Scale Structure and Supernovae la 
combined. 

PACS numbers: 98.80.-k 



Observations of Supernovae la (SNe la) Jlj , the Cos- 
mic Microwave Background (CMB) |2j. [j. |j. [[j] and Large 
Scale Structure (LSS) [E0 

all point towards some form 
of dark energy. Over the years, theorists have come up 
with various models to explain the nature of dark energy 
H H El El El El El (to name a few). In particular, 
it seems plausible that dark energy may be described by 
an (effective El El El El) 

scalar field. Provided no 
symmetry cancels it, there will be a term in the effective 
Lagrangian coupling baryonic matter to the scalar field. 
If this field evolved over cosmological times, such a cou- 
pling would lead to a time dependence of the coupling 
"constants" of baryonic matter El EHH IH El ■ In- 
deed, bounds on the time variation of these "constants" 
restrict the evolution of the scalar field and the strength 
of this coupling [lflj . 

Since the days of Dirac , it has been speculated that 
the fundamental constants of nature may vary. In a real- 
istic GUT scenario, the variation of different couplings is 
interconnected [11 l^l ■ We will ignore this interde- 
pendence in the following and concentrate on a change in 
the fine structure constant a with all other "constants" 
fixed. Lacking detailed knowledge of the dependence of 
a on the scalar field, we make the Ansatz 1 
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Thus, fixing the Taylor coefficient d In a /dip determines 
the evolution of a as a function of p. This has immediate 
consequences: if the fine structure was indeed different at 
higher redshifts, the scalar field must have evolved since 
then. Yet, the Oklo nuclear reactor strongly limits the 
change of a at low redshifts z « 0.1. In our Ansatz, 
freezing a is equivalent to slowing down the evolution 



FIG. 1: Marginalized likelihood distribution for the rel- 
ative change of the fine structure constant quoted as 
log 10 (Aarec/a), so that Aa rcc /a = _io lo sio( At *.W^ The 
distribution is obtained from Set I, i.e. Equivalence Princi- 
ple, Oklo and QSO observations which are all at consider- 
ably lower redshifts than recombination. The dark shaded 
(blue), moderately shaded (green) and light (red) rectan- 
gles correspond to one, two and three a confidence regions. 
The distribution peaks at —4.7 and has width +0.36, —0.26. 
Thus the relative change of a at recombination is Aa ICC /ci = 
-2±2-» x 10" 5 . 



of ip. As the kinetic energy (of a canonical scalar field) 



is given by T = ~a 



this leads to a drop in T and 



therefore to an equation of state w that approaches that 
of a cosmological constant 



T — V 
T + V 



-1. 



(2) 



1 This form has previously been chosen in l2St l29l LiCl ; Clearly, 
an expansio n to higher order or an analytic dependence is also 
conceivable Li4l and may fit the data better. 



In such a scenario, the equation of state of the scalar field 
crosses over from a more positive value at earlier times 
to a behavior that today mimics a cosmological constant 
[In], |H IH ■ I n or< ier to describe such models, we employ 
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FIG. 2: Equation of state for two different parameter sets of 
Equation versus redshift. The equation of state in the 
early Universe is given by w m (here it is set to —0.1), while 
the one today is given by wq (here —0.99). The scale factor 
of the cross-over is determined by a c (0.3 for the solid curve, 
0.7 for the dashed curve), while the width of the transition 
is adjusted by A (0.01 for the solid curve, 0.1 for the dashed 
one). 



the parameterization of pjfjj 

w(a) = W + (w m - Wq) X 



1 + cxp(a c /A) 
1 + cxp([a c — a] /A) 
1 - expOl - a}/ A) 
l-exp(l/A) 



(3) 




FIG. 3: Minus the relative change of the fine structure con- 
stant Aa(z)/ao — [a(z) — ato] /ao vs. redshift quoted as z+l. 
The QSO data is depicted at redshifts z ~ 2 — 3. In addi- 
tion, we indicate the bound from the Oklo natural reactor 
Aa(z m 0.13)/ao > —9 x 10 -8 by a slim vertical box in the 
lower left corner. The solid black curve is a model we picked 
at random from the Monte Carlo and which has parameters 
Q m h 2 = 0.12, h = 0.55, w = -0.95, w m = -0.26, a c = 0.8, 
log 10 (A) = -2.8 and Aa rcc /a = -2.4 x 10" 5 . Please note 
the moderate running of a at early times which is a direct 
result of the scalar field dynamics: with w m = 0, one gets 
ip tx ln(l + z) and hence a moderate evolution at high red- 
shifts. For w m < 0, the running is suppressed further as the 
total energy and hence the available kinetic energy of dark 
energy is considerably lower than for w m = 0. Hence, the 
change of a at nucleosynthesis would have negligible effect 



This versatile parameterization is characterized by the 
equation of state of dark energy today wq, the dark en- 
ergy equation of state during earlier epochs w m , a cross- 
over scale factor a c and a parameter A controlling the 
rapidity of this cross-over (see also Figure |3 for illus- 
tration). It has been shown (and we will confirm this 
later) that some of the four parameters wq, w m ,a c , A 
are not too well constrained by present SNe la, CMB 
and LSS data .'Hi . In principle, a change of a at high 
redshifts will alter recombination. This effect has been 
discussed in and we have modified the REC- 

fast 0] implementation of Cmbeasy 0] to accommo- 
date this deviation from standard recombination. Our 
results agree well with those of [3^, |4^. However, we 
will soon see that Oklo and Quasi Stellar Objects (QSO) 
observations when explained in our framework limit the 
relative change of a at and during recombination. Wc 
will in the following quote this change at the pivotal 
redshift z rcc = 1100. 2 Typically, in our scenario 



2 We identify z = 1100 with the redshift of recombination. This 
is purely semantics and the mismatch with the true redshift of 
recombination is irrelevant. 



Aa rcc /cv = [a(z Tec ) — cto] / ao < 10~ 4 (see also Figure 
|3J|. For such minute changes in a, the recombination 
history remains practically unaltered. Therefore, consid- 
erations like that of |39l l40| have negligible effect. This 
fact allows us to split the analysis in two parts (Set I and 
Set II): Set I with Equivalence principle, Oklo and QSO 
observations for which we do not need to include cosmo- 
logical perturbations. Set II with WMAP, ACBAR, CBI, 
VSA, SDSS and SNe la 0,0 HQ, HQ observations (and 
hence cosmological perturbations), but constant a = ag. 
In both cases, we employed the Markov Chain Monte 
Carlo (MCMC}_p4| from the AnalyzeThis package 
of Cmbeasy |42|. 

For completeness, let us describe the well known recon- 
struction of the scalar field evolution from some parame- 
terization w(a): the starting point is a reconstruction of 
Pde(a) by integrating dlnp^e = — 3[1 + w(a)]dlna. This 
fixes ip (a), because T — |(1 + w)pd e = \a~ 2 (p 2 . From 
this, tp(a) simply follows by integrating ip, where we use 
an integration constant to fix tpo = for convenience. De- 
manding that a should change by Aa rec /a at z lec = 1100 



3 



and using ipo = 0, we can re- write Equation as 

/ \ ^ Aa rcc ^ <p(z) 

a{z) — a + a 



For QSO observations, we use the sample of |46|, |5 



a J tp 



(4) 



For Set I, the parameters of the model were 
il m h 2 , h, w , w m , a c , log 10 (A) and log 10 (Aa ICC /a). 3 
Please note that we omit the spectral index n, optical 
depth r and baryon fraction Slbh 2 which are present in 
the Set Il-analysis, because Set I data does not contain 
information on these parameters. Thus, the distribution 
in n, r and f2t>/i 2 is flat and our simulation for Set I can be 
thought of as trivially marginalized over n, r and fli,h 2 . 
The Equivalence Principle is tested via the differential 
acceleration r\ between test particles of equal mass but 
different composition [2j| 



t)=-3x 1.75- 10" 



d In a 

dip 



AR(1 + Q). (5) 



Here, AR = AZ/(Z + N) « 0.1 is the difference in 
composition and the factor (1 + Q) encodes the theo- 
retical uncertainty in calculating rj. As the true under- 
lying theory is unknown, we 'guesstimate' conservatively 
(1 + Q) € [10~ 2 , 10 2 ] which is used as a flat prior. As 
current experiments give a null result |47| 



ry = ± 3 x 10 



-13 



(6) 



(1 + Q) = 10 -2 (i.e. the lower bound) always agrees best 
with observations. 

Roughly 1.8 billion years ago (i.e. z w 0.1), the Oklo 
natural reactor was up and running. Until recently, the 
reactor was seen as providing yet another null result, 
namely 



Aa 



a oklo 



Otg 2 x 10" 



(7) 



Yet, the latest calculation 50] yields a positive detec- 
tion of 



^ =45+I 5 xlO- 8 
a oklo 



(8) 



To complicate matters further |49j, Fujii et. al [5l| find 
two values, 



Aa 

a oklo 
Aa 

a oklo 



-0.8 ± 1.0) x 10" 
.8 ±0.7) x 10~ 8 . 



(9) 
(10) 



We will take a conservative stand here (given the scat- 
ter in the calculations) and use the null-result of [lif . 
Equation Q. 



which roughly translates into 



Aa 



— = (-5± 1) x 10~ 6 . 
a qso 



(11) 



There have been contradicting claims about such QSO 
results . While time may tell what a really was at 
z w 2, we stick with Equation i|ll|l for the purpose of 
this letter. So we ask: if a really changed in the past, 
what does it tell us about dark energy? 

The main results of our Set I analysis are summarized 
in Figure ^and the left column of Figure |SJ From Figure 
n w e see that Aa rcc /a = —2^ 3 x 10~ 5 . Likewise, in the 
left column of Figure [S] we see that these measurements 
lead to strong bounds on the equation of state today. 
At 2cr confidence level, wo < —0.97. In addition, the 
transition is restricted to occur rather swiftly (as seen 
from log 10 (A)) and at rather recent times (as seen from 
a c ). 

Turning to our Set II analysis, the parameters were the 
amount of matter and baryons fl m h 2 and flbh 2 , Hubble 
constant h, optical depth to re-ionization t and spec- 
tral index n plus four parameters wq, w m , a c , log 10 (A) 
of the dark energy model. As said, we did not need to 
take a varying a into account, because the allowed val- 
ues of Aa rcc /a from Set I are minute in our framework 
and will not alter standard recombination. 5 To test this 
model, we compared to WMAP, VSA ACBAR, CBI, 
SDSS and SNe la measurements Q, S H 0, H @ . The 
versatile (and therefore difficult to constrain) parameter- 
ization of Equation |JSJ has been under investigation in 
[3^|. Our results agree well with those presented in that 
paper (c omp are the right column of our Figure [5] to Fig- 
ure 6 of [36|). With the exception of Wq, the constraints 
are considerably less tight than those inferred from the 
fine structure data. Combining Set I with Set II leads to 
the likelihood distributions shown in Figure 0] and rep- 
resent our main results. At 2a confidence level, we get 
wo < —0.97 and A < 0.1. Quoting the la error bars, we 
get w m = -0.2t° a f Q and a c = 0.79t° o f 2 . 

In this letter, we presented a quantitative analysis of a 
scenario in which the running of the fine structure con- 
stant is driven by a scalar dark energy field. We param- 
eterized the evolution of the scalar field in a versatile 
manner thus covering nearly all of today's quintessence 
models with w > — 1. Together with the Ansatz that 
the induced running of a is linear in the field we found 
stringent constraints for dark energy model building. For 
future work, the analysis may be extended in two direc- 



We picked negative Aa/a, i.e. A« rcc /o = -10 lo sio( A< *«=<=/a) , 



4 To prevent misunderstandings: we used the binned data shown 
in Figure 131 

5 We did perform a Monte Carlo including varying a. As expected, 
nothing changed. 
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FIG. 4: Marginalized likelihood distribution from the combination of Set I and Set II. The dark shaded (blue), moderately 
shaded (green) and light (red) rectangles correspond to one, two and three a confidence regions. At 2<r confidence level, we get 



too < —0.97 and A < 0.1, while with la error bars, we have w m = — 0.2+Q'gQ and a c — 0.79 



+0.06 
0.12- 



tions: as the running of couplings in realistic GUT sce- 
narios is interdependent, one may investigate a running 
of a together with a running of the Planck mass. Sec- 
ondly, the Ansatz of a linear dependence of a on the 
dark energy field may be promoted to higher order or an 
analytic function. In this substantial change of 

recombination or nucleosynthesis physics is well conceiv- 
able. 
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FIG. 5: Marginalized likelihood distribution for the parameters of Equation 0. The left column derives from Oklo, QSO and 
Equivalence Principle observations, while the right column uses WMAP, ACBAR, CBI, VSA, SDSS and SNe la data. Please 
note that the bound on all parameters on the left is considerably tighter than that in the right column. The dark shaded (blue), 
moderately shaded (green) and light (red) rectangles correspond to one, two and three a confidence regions. 



